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The spin- 1/2 quantum Ising chain in a transverse random 
magnetic field is studied by means of the density-matrix renor- 
malization group. The system evolves from an ordered to a 
paramagnetic state as the amplitude of the random field is 
increased. The dependence of the magnetization on a uni- 
form magnetic field in the z direction and the spontaneous 
magnetization as a function of the amplitude of the trans- 
verse random magnetic field are determined. The behavior of 
the spin-spin correlation function both above and at critical- 
ity is studied. The scaling laws for magnetization and cor- 
relation functions are tested against previous numerical and 
renormalization-group results. 

PACS numbers: 75.40.Mg, 75.10.Jm, 05.50.-hq 



I. INTRODUCTION 

The spin- 1/2 quantum Ising model in a transverse ran- 
dom magnetic field (QIMRF) has been the subject of 
an increasing interest in recent times. The main prop- 
erty of the QIMRF is that it undergoes a disorder-driven 
quantum phase transition in the ground state. More- 
over, the critical point of the QIMRF is not an isolated 
point of nonanalyticity, and the zero-field susceptibihty 
diverges in a wide region around the critical point. This 
behavior is due to the existence of diluted domains of 
the wrong phase in the system. Such domains are sta- 
bilized by the disorder. This nonanalytical behavior is 
much more pronounced in the quantum model than in 
its classical counterpart, pointing out the difference be- 
tween quantum and thermal fluctuations in disordered 
systems. More recently the QIMRF has been obtained 
as an effecjtave model for the one-dimensional Kondo lat- 
tice modelEJ (KLM) by means of bosonization of the con- 
duction electrons. The critical behavior of the KLM has 
indeed features which could be explained in terms of the 
corresponding QIMRF .H 

Previous numerical resultaj exploited the mapping of 
the QIMRF in zero external uniform field onto a model 
of noninteracting fermions. This mapping allows for a 
detailed study of correlation functions, but is limited 
by the maximum size of the system that can be han- 
dled. On the other hand, in the presence of an external 
uniform field in the z direction the model is not solv- 
able. In this paper we make use of the density-matrix 
renormalization-group (DMRG) approach in its infinite- 
size version to overcome the limitations imposed by the 



mapping to free fermions. By explicitly introducing an 
external uniform field we are able to study the magne- 
tization and extrapolate the spontaneous magnetization, 
thus testing the scaling results and determining the cor- 
rections beyond scaling. Moreover our computing time 
depends but linearly on the size of the system, for a given 
accuracy. 

In this paper we study the spin-1/2 QIMRF defined by 
the Hamiltonian 



n 



(1) 



where the Sf = o'f/2, {erf} are Pauli spin matrices, 
hi are random-field values uniformly distributed in the 
interval —ho < hi < ho, and hz is a constant field in 
the z direction, added to produce spontaneous symme- 
try breaking in our subsequent calculations. We take the 
units of energy such that J = 1. For the sake of defi- 
niteness we study the ferromagnetic case J > 0, to which 
the antiferromagnetic model can be reduced by a gauge 
transformation on every second site S'f = (— l)*^^ , i be- 
ing the site index. Obviously in the antiferromagnetic 
case hz has to be taken as a staggered field. Moreover, 
a gauge transformation hi — sgn{hi)hi, Sf — sgn{hi)Sf 
can be performed to make all hi positive. 



II. PREVIOUSLY KNOWN ANALYTICAL 
RESULTS 



The QIMRF was studied, analyticallypby McCoy and 
Wu,a Shankar and Murthy,El and Fisher Fisher, in par- 
ticular, makes use of the renormalization group, to deter- 
mine the critical properties of the model defined by the 
Hamiltonian 



J 



^hf 



(2) 



with all jf and hf positive random variables, 
drawn independently from two distributions of densities 
^{J),p{h). 

The phase transition is found to occur at the point 
where \nh^ = InJ^, where the bar denotes average over 
disorder. The distance from criticality is defined as 



5. 



Inh^ ~ In 



var(lnft,^) + var(lnj-'^) ' 



(3) 



In Fisher's notation, the distributions chosen in this 
paper are 
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7r(J) = 5(J-1), pih) 



hg^ for < /i < ho, 



otherwise. 



(4) 



The relations between our and Fisher's variables are J — 
4jf, h, = 2hf, = 2H. 

Both in the ordered iSn< 0) and in the disordered 
{S > 0) phase a GriffithsB region does exist in which 
the magnetization is nonanalytical at iJ = 0. The 
weakly ordered Griffiths "phase" extends over the region 
where mm{jf} < niax{/ij^}, S < 0. The weakly dis- 
ordered Griffiths "phase" extends over the region where 
max{jf } > mm{hf}, S > 0. 

Fisher has found the exact critical scaling function of 
the magnetization and he showed in particular that the 
magnetization at the critical point behaves as 



M {S = 0,H)^ [In (Dh/H) 



(5) 



where = (1 + x/5)/2 = 1.618033... is the golden 
mean and is some nonuniversal scale factor. The 
spontaneous magnetization in the ordered phase depends 
on the distance from criticality according to 



M((5 < 0,i? = 0) 



(6) 



for 6^0 . Further away from the critical point, in the 
weakly ordered phase ((5 < 0), for small H 



M (H) - Mo - [In iD/H)Y 



(7) 



where Mq is the spontaneous magnetization, the ex- 
ponent X cannot be determined within renormalization 
group, D is a nonuniversal scale factor, and 



1 

In the weakly disordered phase ((5 > 0), instead 



MiH) 



{Hy'^[\n{D'/H)] 



1 + 1/2 



(8) 



(9) 



where D' is another nonuniversal scale factor, and H <C 
D'. 

The two-point spin-spin correlation function 



Civ) = {a^oO 
has the typical behavior 

- InC(r) - V? 



(10) 



(11) 



at the critical point (5 = in the limit of large r. The 
(• • •) denotes the quantum average in the ground state. 
In the disordered phase 5 > 0, C(r) should obey the 
scaling equation 



C(r;(5 = 0) 



In 



lnC(r/C) « -r/C, 



(12) 



for r ^ ^, where the typical correlation length ^ has the 
behavior 



(13) 



with V = 1. Young and RiegerB have, by mapping the 
model in Eq. (|l|) for = Q onto a model of free fermions, 
numerically calculated the typical correlation functions 
on finite-size chains. They obtain that their data scale 
best for v = 1.1. 



III. DMRG ALGORITHM 

The DMRG algorithm introduced by Whitei is used 
with minor changes. The random field breaks the sym- 
metry of the chain under reflection and separate density 
matrices are needed for the left and the right subsystems. 
At each step of renormalization a new block is defined by 
truncating the Hilbert space of the left and right sub- 
systems onto the eigenstates of the density matrix of the 
subsystem which correspond to the largest eigenvalues. 
In our calculation the minimum number of states kept is 
8 at each iteration step. However, if a degeneracy occurs 
in the density matrix, then all the degenerate eigenstates 
are kept, a procedure which significantly decreases the 
truncation errorl. Even at the point of the phase tran- 
sition we found a relative truncation error of the order 
10^^*^, which makes our results very accurate. 

The DMRG algorithm is used to calculate the average 
magnetization of the chain and the spin-spin correlation 
functions at some given value of the amplitude of the ran- 
dom field. In the first case the average of the z compo- 
nent of the spin operator of the full system is calculated 
after 50 DMRG steps (the system is then composed of 
102 sites). The value is averaged over a large number 
(usually 1000, sometimes up to 6000, depending on the 
values of hz) of different random-field configurations. In 
the second case = 20 — 40 DMRG steps are ignored to 
make the system large enough so that the edges would 
have a small influence on the innermost spins. One spin 
in the middle of the lattice is marked a and the matrix 
of its z component is at the step iV + 1. The chain 
can be represented as BlgsBi^, where denotes the 
left block of spins, Bj^ the right block and s is another 
spin site in the middle of the lattice. The correlation 
function of the a spin with itself is trivially calculated as 
C(0) = {6l (S'a)^ (E)Ss(g) Sr) = 1/4. The 5's are unit 
matrices of the corresponding parts of the system. At the 
step A^-|-2 the left subsystem B^a is renormalized defin- 
ing a new block B^, according to White's procedurep in 
our procedure we cannot, however, exploit the symmetry 
of the system under reflection, so that the right subsys- 
tem has to be renormalized separately, yielding a new 
block B'j^. Finally two more sites are inserted between 
the new blocks. The system looks now like B'j^bsB'j^, 
where the left one of the two new sites is called b. The 
spin a is now the rightmost spin in the block , so that 
its spin matrix is no longer equal to the simple 2x2 spin 
matrix, but it is some matrix ^^(l) of dimension to x to. 
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obtained by truncating the product Sl ^E) S^. Here m is 
the number of states retained in the block at each step 
(8 or more). The number 1 in the parentheses denotes 
the distance of the spin a from the right edge of the B^. 
The correlation function between the spins a and b is 
C(l) = (Slil) ^ d) 5s^ S'b). After each RG step 
the spin a "moves" deeper into the block Bl and a cor- 
relation function between spins that are farther apart is 
calculated according to C(r) = {Sl{r)iS)S^'SiSs'dS'g). We 
repeat the procedure until a sufficient number (50) of val- 
ues of C(r) is obtained. All of them are averaged over a 
large number of random field configurations (3000-8000). 

It should be mentioned that the magnetization must 
be calculated at nonzero (though very small) to have 
spontaneous symmetry breaking. This small additional 
field makes the magnetization larger than the sponta- 
neous magnetization Mq — M{hz = 0). Several runs at 
different hz are needed to extrapolate Mq. We shall re- 
turn to this point later. On the other hand, as usual, the 
correlation functions are directly calculated at — 0. 

IV. NUMERICAL RESULTS AND DISCUSSION 

We carefully tested the reliability and accuracy of our 
algorithm in the two exactly solvable limits of Hamil- 
tonian (|l|), i.e., Hq = (pure Ising model) and J — 
(noninteracting spins in a random field). In the first 
case we reproduced the exact ground-state wavefunction 
1$) = J]. I (for hz 0+), the ground-state en- 
ergy Eq = —J{N — 1), where N is the number of sites, 
and the correlation function C(r) = 1/4, Vr. In the 
second case we checked that for a given random-field 
configuration, the ground-state wavefunction is 1$) = 
2-N/2 ji _|_ sgn{hi)\ li)] with a ground-state energy 
Ef) = —l/2J2i l^il- We also reproduced the correct aver- 
age of the ground-state energy and correlation functions 
over disorder. 

0.5 I ' , ' , ' , ' , ' 1 




0.1 I ' ' ' ' ' ' ' ' ' 1 

20 40 60 80 100 

RG steps 



FIG. 1. The magnetization dependence on the system 
size. All our further calculations are interrupted after 50 RG 
steps, when the magnetization has saturated. 
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FIG. 2. The dependence of the magnetization M 
and the estimated error of the mean cr^ on the num- 
ber of random magnetic field configurations N^. Here 

(7^ = \jYji{^i - / \J Nc{Nc - 1), where a bar denotes 
the average over configurations. 

The dependence of the calculated quantities both on 
system size (i.e., on the number of RG steps) and on the 
number of random-field configurations has to be carefully 
discussed. A sufficiently large number of RG steps should 
be performed to eliminate the effects of the edges of the 
lattice. As is seen in Fig. |l|, the magnetization is an 
increasing function of the size of the chain. The effect in 
this case becomes typically negligible after 50 RG steps. 
Figure || shows that oscillations of the magnetization as 
well as the estimated error of the mean decrease quite 
rapidly with increasing number of field configurations. 
Usually 1000 configurations are enough for the relative 
error to become 1 — 2 % though a larger number of them 
is required at small external fields. Again we point out 
that in all our calculations the relative truncation error 
was of the order 10~^°, so that the errors that affect 
our data are purely statistical, the convergence to the 
infinite-size limit being controlled by the number of RG 
steps, and the accuracy of our averages being determined 
by the number of random-field configurations. 

The magnetization dependence on the amplitude of the 
random magnetic field /iq at /iz = 1 x 10^^ is shown in 
Fig. ^. The system is in a ferromagnetic state at small 
fields and M = 1/2. The magnetization decreases with 
increasing /iq and the system is driven towards a param- 
agnetic state. Equation (^ predicts the phase transition 
to occur at = e/2 = 1.3591 . . . since the distribution 
(^) gives in our case 

5 = ln^, (14) 

e 

once the mapping to Fisher's variables is taken into ac- 
count. The critical point is fairly well reproduced by our 
extrapolation. 
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FIG. 3. The magnetization dependence on the amph- 
tude of the random magnetic field at the small ordering field 
H — = 5 X 10~® and H = 0. The error bars are smaller 
than the symbols. The data are averaged over 1000 samples. 
The solid line is a fit according to Eq. (^). 

The spontaneous magnetization M{hz = 0) decreases 
according to Eq. (o) . The solid line is the best fit of the 
form Mo = a(-(5)P with a = 0.56 ± 0.02 and (3 = 0.40 ± 
0.02. The exponent (3 agrees well with the predicted value 
2- (j) = (3 - \/5)/2 « 0.382. The values of A/o where 
obtained by fitting Eq. (^. In our fit we find that the 
exponent x is negligibly small. 

The magnetization as a function of the external uni- 
form field at the criticalpoint is plotted in Fig. |^. The 
result agrees with Eq. (||) . The solid line is a fit of the 
form M{H) = c [In {Dh/H)f~^ with c = 0.34 ± 0.01 and 
Dh = 0.024 ± 0.005. 

We also addressed the issue of determining the devia- 
tions of the exponent z from its scaling value 1/(2|(5|). By 
fitting Eq. (0) and Eq. (|) for <5 < and J > 0, respec- 
tively, we obtained the 1/z as a function of 2\S\, plotted 
in Fig. |5[ The exponent z is smaller than its scaling value 
in the ordered phase, and larger in the disordered phase. 
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FIG. 4. The magnetization dependence on the ordering 
field at criticality. The solid line is a fit of the form (^. 

The average of the logarithm of the two-point spin- 
spin correlation function at the critical point is shown in 
Fig. 0. According to Eq. ([ll|) , InC should be proportional 
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FIG. 5. The dependence of the inverse of the exponent z 
on the distance from criticality. The solid line is the theoret- 
ical behavior oi 1/z for small 5 values, according to Eq. (^). 

to the square root of r for large r. The solid straight line 
is a fit to the data with 20 < r < 50 and has a slope 
-0.96 ± 0.01 and an intercept 0.52 ± 0.01. Our result is 
thus 



InC = 0.96(1) 



.1/2 



0.52(1), 



(15) 



where the figures in the parentheses denote the standard 
error. 
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FIG. 6. The average of the logarithm of the correlation 
function at the critical point. The straight solid line is a fit 
to the data with 20 < r < 50, the slope of 0.96 ± 0.01, and 
the intercept of 0.52 ± 0.01. The logarithm of the correlation 
function was averaged over 8000 random field configurations. 

Our best scaling plot of the logarithm of the correlation 
function in the disordered phase is shown in Fig. 0. The 
data for different 6 collapse on a straight line, of slope 
« —1, when the exponent is iy = 1.04 which is slightly 
larger than the scaling value, predicted by Fisherjj v — 
1, but smaller than the exponent determined by Young 
and Rieger,H which obtained the best scaling at D = 1.1. 
That the slope of the line is « — 1 lends strong support 
to the suggestion by Fisher, mentioned in Ref. ||, that 
= \nh^ — InJ^ , which for the distribution (4) gives 
that = d 01 ln(7(r/f) « -6r for r > |. 
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FIG. 7. A scaling plot of the average of the logarithm 
of the correlation function in the disordered phase with the 
exponent v = 1.04. 



V. CONCLUSIONS 



optimal doping could be determined by the physics ofj-a 
nearby quantum critical point, associated to a magnetic L£ 
or an incommensurate charge-density-wave instability.ti 
The anomalies observed in the optimally doped metallic 
phase and the behavior of the critical temperature have a 
natural interpretation once the lack of any characteristic 
energy scale beside the temperature itself is assumed as a 
consequence of the presence of a quantum critical point. 

The approach described in this paper allows for a 
straightforward application to all those quantum systems 
that can be mapped onto effective spin models in the 
presence of random fields. Our results show indeed that 
the DMRG algorithm is a good and reliable alternative 
to other numerical methods within two respects: the pos- 
sibility of studying systems of a considerable size, when 
needed, and the possibility of introducing an external 
field coupled to the order parameter, to extrapolate its 
behavior in the state of spontaneously broken symmetry. 



In summary, we have analyzed some properties of the 
spin-1/2 quantum Ising model in a transverse random 
magnetic field (QIMRF) by means of the density-matrix 
renormalization group (DMRG). This method has two 
advantages with respect to other numerical methods. 
First of all, the computing time, for a given accuracy, 
scales linearly with the size of the system. Thus there is 
no severe limitation in increasing the size of the system. 
Our results show that typically the algorithm converges 
to the infinite-size limit after 50 - 100 DMRG steps, i.e., 
when the system is composed of 100 - 200 sites. Second, 
an external uniform field in the z direction could be ex- 
plicitly added, allowing us to determine the scaling laws 
for the magnetization^ both as a function of the external 
field and of the distance from criticality. This is the most 
relevant aspect of our approach. Indeed the introduction 
of the external field spoils the integrability which was 
exploited in Ref. ||, but is absolutely harmless from the 
point of view of the DMRG, the only effect being that 
a large number of random-field configurations is needed 
for a given accuracy. 

It is evident that once the solvability of the model is not 
required, other systems, which cannot be mapped to free 
fermions, may equally well be investigated by means of 
the DMRG. The relevance of our results stands then not 
only in the accurate test of predicted scaling relations, 
and in the determination of nonuniversal coefficients and 
corrections beyond scaling, but mostly in the wide appli- 
cability of the method described above to the problem of 
quantum phase transitions in disordered systems. This 
is a relevant issue in solid-state physics. For instance, 
as we discussed in the Introduction, the QIMRF was ott 
tained as an effective model for the Kondo lattice model.El 
Moreover, recent theoretical developments point out the 
relevance of quantum phase transitionsEI to the physics 
of high- Tc superconductors. It has indeed been recently 
suggested that the properties of these systems close to 
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